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UHE cosmic neutrino interaction with the cosmic neutrino background (CnuB) is expected to produce ab- 
sorption dips in the UHE neutrino flux at energies above the threshold for Z-boson resonant production. The 
observation of these dips would constitute an evidence for the existence of the CnuB; they could also be used 
to determine the value of the relic neutrino masses as well as some features of the population of UHE neutrino 
sources. After breafly discussing the current prospects for relic neutrino spectroscopy, we present a calculation of 
the UHE neutrino transmission probability based on finite-temperature field theory which takes into account the 
thermal motion of the relic neutrinos. We then compare our results with the approximate expressions existing in 
the literature and discuss the influence of thermal effects on the absorption dips in the context of realistic UHE 
neutrino fluxes and favoured neutrino mass schemes. 



1. Introduction and motivations 

The existence of a cosmological background of 
relic neutrinos (CVB), characterized by a present 
temperature of T„ ~ 1.95 K (1.69 x 1CT 4 eV) and 
a number density n„ m 56 cm -3 per species, is 
a robust prediction of Big Bang cosmology pQ. 
Its presence, at least at early times, seems now 
confirmed by data coming from observational cos- 
mology, which also provide bounds on the num- 
ber of neutrino species (now greater than zero) 
and on their absolute mass scale (see for an 
extensive discussion of these issues). Direct de- 
tection of the CVB at present time is however 
much more problematic. Several laboratory ex- 
periments (discussed i.e. in [3]) have been pro- 
posed to detect relic neutrino through their weak 
interactions, using torsion balances or even ac- 
celerated beams of nuclei, but the technological 
improvements required in those experimental se- 
tups are way beyond the state of the art. 

An interesting alternative is provided by ex- 
ploiting cosmic rays - actually ultra-high energy 
cosmic neutrinos (UHEj/) - as a natural beam and 
search for evidence of their interaction with the 
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relic neutrino background. Two approaches have 
emerged in that context: the first one searches 
for "emission" features and the production of 
charged cosmic rays or photons beyond the GZK 
cutoff, through the so-called "Z-burst" mecha- 
nism 1 113) (the dominant UHE^- Cz^B interaction 
channel being v v — > Z — > X). This mecha- 
nism probes the universe delimited by the GZK 
sphere and is therefore sensitive to local overden- 
sitics in the Ci/B that may arise from the gravi- 
tational clustering on massive structures, like the 
Virgo cluster [S], and would result in a directional 
excess of UHE cosmic rays. 

The second approach, on which we will con- 
centrate here, looks for "absorption" features in 
the UHE^ flux that would reflect their interac- 
tion with the relic neutrinos along their path 
[7181911011 J] . As illustrated below, the position of 
the absorption dip approximately corresponds to 
the rcdshiftcd resonance energy for the Z-boson 
exchange, namely 



so that it directly depends on the value of the neu- 
trino mass m v and of the source redshift z s . The 
shape and depth of the dip also reflect to some 
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extent the neutrino mixing pattern as well as the 
characteristics and distribution of UHE^ sources. 
In that sense, and provided adequate sensitivity 
and energy resolution of the detectors, a detailed 
spectroscopic study of the UHE^ flux could allow 
the determination of the absolute neutrino masses 
along with an experimental proof of existence of 
the Cz^B |12I13| . This mechanism has generated 
a renewed interest in the light of the recent re- 
sults obtained both in observational cosmology 
and for determination of the neutrino mass and 
mixing schemes. 

However, most of the work in the literature de- 
scribe the UHEz/-CVB interactions assuming that 
relic neutrinos are at rest, while the effects of 
thermal motion should be included as soon as 
the momentum of relic neutrinos becomes com- 
parable to their mass (or even before) . Therefore 
we have chosen to address the question of relic 
neutrino spectroscopy in the framework of finite- 
temperature field theory (FTFT), which allows 
to take into account the thermal effects in a sys- 
tematic way {Tj\. In sectional we briefly present 
our calculation of the damping of an UHE neu- 
trino travelling across the Cz^B and determine the 
absorption probability for a neutrino emitted at 
a given redshift. We then illustrate our calcula- 
tions in two realistic contexts and explore various 
combinations of parameters to investigate the dif- 
ferences between the FTFT calculation and previ- 
ous approximations. In section [3] we discuss how 
the thermal broadening of the absorption lines in 
the UHEi/ flux could affect the determination of 
m v and of the characteristics of the UHEz/ source 
population. In section we comment our results 
in the context of relic neutrino clustering, for dif- 
ferent hypothesis on the density and scale of the 
clusters. Concluding remarks are presented in 
section 

2. Damping of UHE neutrinos across the 
relic neutrino background 

The equation of motion of an UHE neutrino 
with four- momentum fc M = {£ K , K) and mass m v 
traveling across the Cz^B reads 

Cfc- m v - E)^ = (2) 
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Fi gure 1. Feynman diagram for the one- loop self-energy 
of an UHE neutrino due to a Z-boson exchange with an 
(anti-)neutrino from the relic background; the blob on the 
Z propagator indicates that we use the dressed propagator 
and the cut is to select the imaginary part of the diagram. 



where the self-energy E accounts for the interac- 
tions with the surrounding medium. The domi- 
nant process is here the Z-boson exchange in the 
s-channel, as shown in fig. ^ We determine E 
from a FTFT one-loop calculation carried out in 
terms of the (vacuum) Z propagator and the ther- 
mal propagator of the relic neutrinos. The latter 
depends on the functions f v (P) and fa{P) which 
describe the momentum distributions of neutri- 
nos (antincutrinos) in the thermal bath. These 
functions take the simple relativistic Fermi-Dirac 
form 

where T v is the temperature of the QvQ and we 
have neglected the chemical potential. It is worth 
noting here that, although relic neutrinos are not 
relativistic anymore at present time, their distri- 
bution has maintained the form it had at the time 
of neutrino decoupling, corresponding to T v ~ 1 
MeV. 

The dispersion relation corresponding to eq. 
(0) is given by £ K = S r (K) - ij(K)/2, where 
£k and 7 are functions of K. The damping fac- 
tor 7 governs the propagation of the UHEi/ across 
the background of relic neutrinos and is directly 
related to the imaginary part of the self-energy, 
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Figure 2. Top : cross-section a v c{P,K), in units of 10 _31 cm 2 , as given by eq. (6), as a function of the energy of the 
incident neutrino, K, and of the relic neutrino momentum, P. From left to right, the panels correspond to a neutrino mass 
10 — 1 , 10 — 2 , and 10 — 3 eV. Bottom : Transmission probability P^{Kq, z s ) as a function of the UHEi/ energy as detected 
on Earth, Kq, for a source located at redshifts z s = 1, 5, 10, 20 (from top to bottom in each panel) and for a neutrino mass 
m v = 10 _1 , 10 -2 , 10~ 3 eV (from left to right). The continued, black curves correspond to the full damping as given by 
eqs. (5) and (6), while the dotted (red) curves are for the approximation of relic neutrinos at rest, eq.(7). 



Ej I n the approximation that the UHEi/ 

are ultrarelativistic and that we can neglect the 
background effects on their energy (£ r (K) ~ K), 
the damping can be written as (see ^3 f° r the 
detailed calculation) 
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For m v <C Mz , K and neglecting terms of order 
r|/M|, we have 
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where E p = \J P 2 + m 2 is the energy of the 
relic neutrino. Taking the limit of eq. 10 for 



P — > 0, one recovers the approximated cross- 
section used for relic neutrinos at rest, with the 
Z peak at the UHE^ "bare" resonance energy 
M z /(2m v ). The damping reads in that 
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in agreement with the results of JS]. The ex- 
pressions used in can be obtained by fur- 
ther evaluating the cross-section at the pole of 
the resonance, 2mK res = M z (narrow-width ap- 
proximation). However, this approximation is no 
longer valid for relatively small m„ and/or large 
GvQ temperatures: fig. 1 (top line) shows how 
the resonance peak in the vv — > Z cross-section 
broadens and shifts to lower UHE^ energies as 
P increases. Actually two effects combine: the 
modification of the cross-section peak due to its 
dependance in E p , and the thermal distribution 
which selects a range of relic neutrino momenta 
close to the temperature of the Czd3. The net ef- 
fect of this thermal broadening is a reduction of 
the damping, which affects the transmission prob- 
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ability and the depth and shape of the absorption 
dips. 

The transmission probability for an UHEz/ 
emitted at a redshift z s to be detected on Earth 
with an energy Kq is obtained by integrating the 
damping along the UHEz/ path, taking into ac- 
count that both the UHEz/ energy and the Cz/B 
temperature are redshifted: 



Pt(K ,z, 
exp 



dz 



o H(z)(l + z) 



(8) 



where H = H ^0.3(1 + ,z) 3 + 0.7 is the Hub- 
ble factor as suggested by recent observations 
|14| . Fig. 1 (bottom line) compares the trans- 
mission probabilities obtained from eqs. (JHJ) and 
(0), for an UHEz/ emitted at different rcdshifts 
and m v ranging from lO -1 to 1CP 3 eV. As long 
as m u /T v > 1CP 2 , the shape of the absorp- 
tion dip is not affected by thermal broadening 
and is rather sharply delimited, at high energies, 
by the bare resonant energy for the propagat- 
ing neutrino, Kq — K res = Af|/(2m„), and at 
low energies by the redshifted resonant energy 
Kq = K res /(l + z s ). Evaluating the position of 
these points would in principle allow us to deter- 
mine m v as well as z s , the redshift at which the 
UHEz/ was emitted. As m v jT v decreases, how- 
ever, the absorption dips get broadened and shift 
to lower energies, and this effect also increases 
with the redshift since UHE neutrinos from dis- 
tant sources are emitted in a hotter background. 

3. Absorption lines in the UHE neutrino 
flux 

The results presented so far deal with a mo- 
noenergetic source of UHEz/ located at a given 
redshift. A realistic situation would more prob- 
ably involve a distribution of sources which emit 
UHEz/ with some given energy spectrum and fla- 
vor composition, the latter evolving along the 
neutrino pathway. To investigate these effects, 
and following the approach of we have con- 
sidered a flux of UHEz/ of the form: 



dz 

— P T (K ,z) V (z) J v (K ) 



The distribution of sources, 
r)(z) = 7? (1 + z) n 8{z- z mi 



l )0(z u 



*), (9) 



is well-suited for an approximate description of 
models ranging from astrophysical production 
sites ("bottom- up" mechanisms for which n ~ 4 
and Zmax < 10) to exotic, non-accelerator sources 
(which typically have n ~ 1 — 2 and may extend to 
larger redshifts). The injection spectrum is taken 
as a power-law with a cutoff at some high energy 
K m ax > K res {\ + z) (we do not consider here the 
possibility of broken power-law spectra): 



J„(K) = j v K~ a 0(K n 



K) 



(10) 
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with the spectral index a ranging between 1 and 
2, depending on the production mechanism con- 
sidered. Under these assumptions, the normal- 
ized flux only depends on the difference n — a of 
the spectral indexes. We then computed numeri- 
cally the UHEz/ flux produced for two distinct val- 
ues of the spectral index combination: n — a = 2, 
which is representative of a distribution of astro- 
physical sources, and n — a — 0, which could 
describe the UHEz/ flux produced through some 
top-down process. 

As for the flavor content of T v , the ratio pre- 
dicted for standard, hadronic sources is J„ c : 
Jv M ■ Ju T = 1:2:0, while more exotic mech- 
anisms could be characterized by a democratic 
J„ c : J v : J„ T =1:1:1. In both cases, and 
in view of our current knowledge of the neutrino 
mixing parameters, it is reasonable to think that 
the total neutrino flux detected at Earth will just 
be the sum on all three neutrino mass eigentstates 
|l()j . Therefore we considered mass patterns com- 
patibles with the favoured 3-neutrino hierarchical 
schemes ^Hl, namely: 

normal : 1113 = 5 10~ 2 eV, rri2 — 9 10~ 3 eV 
inverted : = = 5 10~ 2 eV 

and chose the values ui\ = 10~ 3 eV or 10~ 4 eV 
for the third, unknown mass. 

The results are presented in figs. |3]and 4, which 
display the all-flavour UHEz/ flux as a function 
of the present energy Kq of the UHE neutrino, 
normalized to the flux in absence of absorption 
effects. Both figures show how thermal broaden- 
ing globally modify the shape and extension of the 




Figure 3. UHEf flux in presence of damping, T v , with a — n = and z s = 10, 20 (typical top-down model), normalized 
to the corresponding flux in absence of interactions. The top row is for the normal neutrino mass hierarchy and the bottom 
one for the inverted hierarchy. Plots on the left have m\ = 10 — 3 eV and plots on the right mj = 10 — 4 eV. Colour code is 
as in fig. 1. 
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Figure 4. UHEf flux in presence of damping, T v , with a — n = 2 and z s = 5, 10 (typical bottom-up model), normalized 
to the corresponding flux in absence of interactions. The top row is for the normal neutrino mass hierarchy and the bottom 
one for the inverted hierarchy. Plots on the left have mi = 10 -3 eV and plots on the right mi = 10 — 4 eV. Colour code is 
as in fig. 1. 
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absorption lines respect to the approximated case. 
The dip corresponding to the smallest mass is al- 
most always washed out and only contribute to 
further broadening the line at high energies. For 
very large z max , the merging of the two other dips 
in the normal hierarchy scheme makes it more dif- 
ficult to differentiate from the inverted one. 

4. Absorption lines in the case of relic neu- 
trino clustering 

The possible clustering of relic neutrinos onto 
dark matter halos has been intensively studied, 
in particular in the context of the generation of 
the UHE cosmic rays through the Z-burst mech- 
anism 1 115 j . Recent calculations using Vlasov 
(Boltzann collisionless) equation have presented 
revised estimations of the density profiles and 
typical spatial extension of the neutrino clus- 
ters |18I19| . They give overdensities of the or- 
der of 10-10 4 no„ and cluster scales L ~ 0.01 — 1 
Mpc, depending on the neutrino mass, the mass 
of the attracting halo and its velocity dispersion 
(typically 200 km/s for a galaxy and 1000 
km/s for a galaxy cluster). Limits to the clus- 
tering of neutrinos on large scales are also set by 
the Pauli exclusion principle and the limit on the 
maximum phase-space density jl'Oj . which imply 
that only neutrinos with mass m„ > 1 eV will ef- 
ficiently cluster on galactic halos (Lq ~ 50 kpc), 
while neutrinos with m„ > 0.1 eV can cluster 
on the much bigger scales associated to halos of 
(super-)clusters (L c ~ 1 Mpc) |5I18I19| . 

To compute the UHEi/ absorption due to clus- 
tered neutrinos, we substituted f v (P) in eq. (JHJ) by 
a modified Fermi-Dirac distribution, 

1 P -*IT„ , i 
f cl (P) = - — I (11) 

which parametrises reasonably well the distribu- 
tions functions presented in |19j in function of a 
single parameter $. The neutrino density corre- 
sponding to ea II II is 

n? = -^(l + e-*/^i 3 (-e*/ T "), (12) 
where Li3 (x) is the trilogarithm function. We also 



assume that the cluster density is constant, i.e. 

f n v = N c i n u0 r < L cl 
\ n u = n v0 r > L ch 

so that f v {P) does not depend on the position. 
For a given overdensity factor N c \ we then solve 
n^} = Ncitivo numerically for 

We then computed the transmission probabil- 
ity for a cluster of relic neutrinos with mass 0.1 
or 1 eV, located between the UHEf source and 
the observer. As expected, the effect of the ther- 
mal motion of the neutrinos is generally negligible 
or small due to the relatively small overdensities 
achievable and the absence of redshift effect in T v . 
Fig. ISJshows that we have to saturate the bounds 
on both parameters N cl and L c \ to obtain a sig- 
nificant effect on the depth of the absorption line: 
for a maximal overdensity factor N c \ — 10 4 , the 
maximum absorption probability across the clus- 
ter is reduced from w 55% to w 35%. 

5. Conclusions 

The study of both absorption and emission fea- 
tures in the spectrum of UHE neutrinos offers 
probably today's most promising method for a 
relatively direct observation of the cosmic neu- 
trino background. This is especially true in view 
of the upcoming generation of detectors, like ICE- 
CUBE [ZD, ANITA j22J, FORTE gSj, GLUE 
|24j . the Pierre Auger Observatory [2H|v which 
are already putting limits on the UHEz/ flux and 
starting to constrain the corresponding theoreti- 
cal models. 

One has however to keep in mind that the ther- 
mal motion of the the relic neutrinos could sig- 
nificantly affect the shape and position of the 
absorption dips in the UHEi/ spectrum. From 
the exploration of the parameter space currently 
allowed by astrophysical and cosmological con- 
straints, we have seen that, even in the regime of 
non-relativistic neutrinos, the dips can be broad- 
ened and shifted to lower energies. This will com- 
plicate their observation in real experiments, es- 
pecially if the UHE^ source population is concen- 
trated at small redshifts, producing rather shal- 
low and extended dips. The shift to lower en- 
ergies, where neutrino fluxes are expected to be 
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Figure 5. Transmission probability for a cluster of extension 1 Mpc, made of neutrinos of mass 1 eV (left) and 10 -1 eV 
(right), with a constant neutrino density n^} = 10 3 n$ v and n^} = 10 4 no„ (from top to bottom in each plot). The colour 
code is the same as in fig. 1. 



higher, could in principle increase the detection 
potential respect to the case of a relic neutrino 
at rest with the same mass. Still, the situation 
will be more intricate if the pattern of neutrino 
mass eigenstates is such that their combined ef- 
fect results in a superposition of dips with differ- 
ent depths and extensions. Considering the small 
number of events expected in the UHEi/ detectors 
at those high energies, it seems rather unrealistic 
to resolve the detail of these complex absorption 
patterns. 

As for the case of UHE^ going through relic 
neutrino clusters, we have seen that the effect 
on the absorption is relatively small compared 
to that of UHE^ that travel on cosmological dis- 
tances, since clustering only occurs at small red- 
shifts. For the same reason, and because only 
rather heavy neutrinos do cluster, thermal effects 
have a limited impact on the shape of the absorp- 
tion lines and only result in their attenuation. On 
the other hand, nearby clusters of relic neutrinos 
are expected to play an important role in the con- 
text of the Z-burst mechanism and the emission 
of UHE cosmic rays. 
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